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I. INTRODUCTION 



Dirac optical models are widely used in analyzing e 



(e, e'p) reactions 
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ectron scattering data for (e, e' 
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Recently these models have shown to give excellent agreement with spin observables for 
16 0(e, e'p) (25! where the behavior of the observables at large missing momentum has been 
attributed to dynamical relativistic effects due to "spinor distortion" . These models have also 
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28j for indications of medium 



been used in the analysis of recent data for 4 He(e, e'p) (26, 
modification of nucleons in the nuclei. Evidence for such modifications in this case relies 
on use of optical model calculations with and without medium modified from factors. Since 
the size of the difference between the calculation without medium modified form factors and 
the data is on the order of 5 to 10 percent, any conclusion based on this approach requires 
that the model calculations can be trusted to a similar level of accuracy. (It should be noted 
that this later effect has been described in a more traditional approach 2^| by including 
charge exchange interactions. The choice of the parameters for are reasonable, but are not 
well constrained by data.) 

Given the importance of the questions that are being addressed with these models, it is 
necessary to consider their foundations. The fundamental assumption is that the nucleon- 
nucleus interaction can be described in terms of a single-particle hamiltonian of the form 



H = -a ■ V + /3 [m + S p T (r, E)\ + V OPT (r, E) 



(1) 



where Sopt(^,E) and Vopt(t,E) are complex, energy-dependent, scalar and vector opti- 
cal potentials. This approach was first used to provide a phenomenological description of 



proton-nucleus elastic scattering 
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31 



, 3^| . It was subsequently demonstrated that optical 
potentials derived from parameterized NN interactions in the impulse approximation pro- 
vided a very good description of the spin observables for proton-nucleus elastic scattering 
at intermediate energies with a minimal number of parameters^. H . I35 ] . It should be 
noted, however, that although the origins of the optical potential in nonrelativistic multiple- 
scattering theory have received a great deal of theoretical attention, the Dirac optical model 
proceeded by analogy to the nonrelativistic case without reference to a relativistic many- 
body theory. 

Similarly, the first applications of the Dirac optical model to (e, e'p) and (e, e') reactions Q, 
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assumed that the necessary current matrix elements could be obtained in analogy 
to the nonrelativistic case with wave functions obtained from one-body Dirac equations 
and the current operator described by a one-body current. In both the nonrelativistic 
and relativistic cases this assumption leads to a lack of current conservation. This lack of 
current conservation is a direct result of the underlying many-body nature of these reactions. 
This manifests itself in several related ways. The first is associated with the composite 
nature of the nucleon resulting in the need for form factors which interfere with the usual 
single-particle Ward-Takahashi identities and imply that one-body current be of a much 
more complicated general off-shell form. The second is associated the with the appearance 
of many-body exchange or interaction currents. Finally, the use of an optical potential 
implies that the many-body problem has been reduced to an effective two-body where the 
contributions of channels associated with excitation of the residual system are hidden in the 
optical potentials. Since these hidden channels can be excited by virtual photon absorption, 
a consistent treatment of the reaction requires that an effective current operator be used in 
place of the simple one-body current. The first source of current non- conservation has been 
addressed by studying the effect of various on-shell equivalent forms of the single-nucleon 
current on the optical model calculations Q, 0, lid ]. This gives some rough indication 
of the size of violation of current conservation, but does not really address the underlying 
problem. The second source has been addressed by including two-body meson-exchange 
currents in an approximate fashion j^]. The problem of the effect of the reduction of the 
many-body problem to an effective optical model on the current has been discussed in a 
general fashion, but has not been studied in any concrete way. jj, [jf| 

The purpose of this paper is to show that it is indeed possible to obtain a Dirac optical 
model from an underlying covariant theory and to obtain the corresponding effective current 
operator necessary to maintain electromagnetic current conservation. In doing this several 
choices will be made in the reorganization of the covariant theory into the optical model. 
Clearly, this approach is not necessarily unique. Therefore, the hope is that this work 
will stimulate the development of alternate approaches with the hope that this will lead 
to an improvement in the phenomenology for the application of the Dirac optical model to 
electromagnetic processes. 

The starting point for this work are the many-body Bethe-Salpeter equations. These 
equations are most easily understood as a resummation of all Feynman diagrams for n-point 
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functions. The n/2 particles associated with the external legs of the n-point function are 
treated as explicit degrees of freedom while all other degrees of freedom are collected into 
a set of irreducible kernels. These degrees of freedom are implicit. The kernels are then 
used in integral equations to sum all contributions to the n-point functions. Since these 
equations are based in Feynman perturbation theory, all elements of the integral equations 
are manifestly covariant. For spin-1/2 constituents the one-body propagators appearing in 
the integral equations are solutions to the Dirac equation so it reasonable to believe that it is 
be possible to reduce the many-body problem to an effective theory involving the interaction 
of a Dirac particle with an (n — l)-body system. The structure of the integral equation for 
the Bethe-Salpeter n-point functions is similar in form to those for nonrelativistic multiple- 
scattering theory with the exceptions that all integrals are four-dimensional rather than 
three dimensional, and that the all propagators are local, whereas propagators of time- 
ordered description usually used in multiple-scattering theory are global. 

For simplicity, the simplest illustrative case of the process, the three-body Bethe-Salpeter 
equation (221, for distinguishable particles, is used to show how the reduction to an ef- 
fective optical model can be implemented. This equation is relatively simple in structure 
and the construction of electromagnetic current matrix elements for this equation is well 



understood 
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411 ]. The optical model is obtained by reducing the three-body problem 



to an effective two-body problem. The effective kernel for the interaction between the bound 
state of two of the particles and the remaining particle can then be interpreted as an "optical 
potential." A similar reduction of the Bethe-Salpeter current matrix elements leads to the 
identification of an effective current operator consistent with the optical model. This optical 
model current will be shown to result in conserved current matrix elements. 

In the first section the optical model for the interaction of one particle with a bound state 
of the remaining pair is constructed. Next, bound and scattering states are defined in terms 
of the optical model states. Finally, the effective optical model current is constructed and the 
impulse approximation contribution to the effective current is isolated. It is then shown that 
the optical model current satisfies a Ward-Takahashi identity involving the optical potential 
which results in conserved current matrix elements. 
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II. OPTICAL MODEL REPRESENTATION OF THE THREE-BODY SCATTER- 
ING MATRIX 



Here we will use a matrix form for the three-body Bethe-Salpeter equation described in 
4l| to simplify the reduction of the three-body problem to the effective two-body problem. 
This is formulation is summarized in the appendix for the convenience of the reader. For 
three distinguishable particles, the three-body scattering matrix can be written in matrix 
form as using (jA5|) 

T = M- MG° BS BT, (2) 

where the matrices are defined in the appendix. 

Our objective is to reduce this expression so than we can extract an effective equation 
for particle 1 scattering from a bound state of particles 2 and 3. This is accomplished by 
separating the two-body t matrix for particles 2 and 3 into terms containing bound-state 
poles and a residual piece containing only the scattering cuts. Assuming that there is only a 
single bound state for particles 2 and 3, the two-body scattering matrix in momentum space 
has the form 

i / r< 2 ) V, i^r^V, p 1 ) -p 1 ', -p^r^i-p 1 , -P l ) 



2E(P) \ pi° - EiP^+irj P l0 + EiP 1 ) -it] 

+M^(p 1 ',p\P 1 ), 



where P 1 is the total four-momentum of the pair, p l and p 1 ' are the initial and final relative 
four-momenta of the pair, T^ 1 is the bound-state vertex function for the pair, m 1 is the 
mass of the bound state and M* is the residual scattering matrix. For relativistic many-body 
equations there is no clean factorization of the vertex functions into relative and center-of- 
mass pieces. This means that vertex functions are explicitly dependent upon the total 
momentum. Any spinor indices associate with the vertex function are suppressed and are 
assumed to be summed. Note that there are positive and negative poles associated with the 
positive- and negative-energy bound states. There are several ways to precede at this point 
in reducing the three-body problem to an effective two-body problem. Both the positive 
and negative poles can be retained and thus explicitly include the interaction with particle 
1 and the negative-energy bound state. This has the virtue that the decomposition of the 
scattering matrix can be written in a manifestly covariant form. However, any attempt to 
write this in the form of an optical model will result in a form which is more complicated 



than is usually assumed. In addition, it is reasonable to assume that the contributions from 
the negative-energy pole will be small especially when this approach is extended to systems 
with more particles which means that the additional complexity may have little real physical 
impact. 

A second approach would then be to treat only the positive-energy pole with the negative 
energy pole becoming part of the residual scattering matrix. This decomposition will not 
be manifestly covariant. Furthermore, the vertex functions are only uniquely defined at the 
pole and (J3J) assumes that they are defined at this point as is indicated by the hat over the 
total momentum. The momenta in (j2J are not similarly restricted and the action of inverse 
two-body propagators on the bound state are not generally defined. 

A third alternative, which is the one used here, is to assume that in any loop integrals 
involving the bound state the positive energy pole will be picked up resulting which will 
restrict P 1 to be on-shell. This prescription is manifestly covariant and the resulting de- 
composition of the two-body scattering matrix is also covariant. This can be realized in the 
equations for the three-body scattering matrix by writing 

M = -D 1 |r^ x ) iG^iiQ 1 ) (r^ 1 ! D lT + Mr = M p + M r , (4) 

where 



D 1 



(5) 
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and Q 1 is an operator that places the total momentum for particles 2 and 3 on the bound- 
state mass shell by requiring that the appropriate pole be picked up in any loops containing 
the two-body t-matrix for particles 2 and 3. 

We can now decompose the equation for the t-matrix according to whether the initial 
(final) interaction is the pole contribution (superscript P) or the residual contribution (su- 
perscript R.) This leads to the set of coupled matrix equations 

T pp = M P - MpG° bs B (T pp + T RP ) (6) 

r RP = -m r g° bs b (r pp + r RP ) (7) 
r PR = -MpG° bs b (r PR + r RR ) (8) 

T RR = M R - MrG° bs B (T pr + T RR ) . (9) 



The solution of this set of equations is facilitated by defining a scattering matrix that does 
not contain any contributions from the bound state for particles 2 and 3. This is defined as 

T R = M R - MrG bs BT r = M r - TrG° bs BMr . (10) 

The second of these two forms can be solved to give 

Mr = Tr + TrG° bs BMr . (11) 

Using this, (JTj) can be rewritten as 

r RP = -(t r + t r g bs bMr) g bs b (r pp + r RP ) 

= -t r g bs b (t pp + r RP ) - TrG bs bm r g bs b (t pp + r RP ) 

= -T R G BS B (T pp + T RP ) + T R G BS BT RP 

= -T R G BS BT pp } (12) 

where we have used the fact that D lT BD 1 = in simplifying the equations. This can be 
used in (jHJ) to give 

T pp = M P - M P G BS BT RP = M P + M p G bs BTrG bs BT pp . (13) 

From (gj), 

M P = —D 1 (r^ 1 ) zGTW) <r (2)1 | D lT = -iG-^D 1 \Y^) (-iG^iQ 1 (T^ 1 \ D lT iG^ . 

(14) 

Using this and iterating of (fTB^) it is possible to make the conjecture that 

T pp = -iG^D 1 ^^iQ't-iGoPTW 1 {T^\D lT iG^ . (15) 

Substituting this into (|T3|) . 

-tG^D^T^iQ^-iGoPTWiT^D^iG^ 1 
= -iG^D 1 \Y^) {-iG^iQ 1 {T^\D lT iG^ 

+iG^D l \T® X ) {-iG^iQ 1 {T^D^iG^GlsBTpGlsBiG^D 1 \T^)iQ l 
x {-iGopt)iQ 1 <r (2)1 1 D lT iG^ 1 . (16) 
We can then identify 

iQ l G PTiQ l = GnQ 1 

-GiiQ 1 (T^\ D^iG^GIsBTrGIsBiG-^D 1 \Y® 1 ) iQ 1 G OP tiQ 1 
= GnQ 1 - GnQ 1 <r (2)1 1 D^GbsBTrGbsBD 1 |r (2)1 ) Q'GopM 1 (17) 



where G BS = — iG^G^ is the free two-body propagator for particles 2 and 3, and l^ 2 ) 1 ^) = 
G BS Jr^ 1 ) is the bound-state Bethe-Salpeter wave function for particles 2 and 3. Defining 

Vqpt = 1Q 1 <$ (2)1 | D lT BT R BD l Q 1 (18) 

as the optical potential, the optical model propagator can be rewritten as 

iQ}G OP TiQ} = (Gi - G^qptGopt) iQ} ■ (19) 

Note that keeping only the leading terms in [TBI yields 

Vop T = tQ 1 (& 2)1 \ {M° + M 2 + M 3 ) I^Q 1 . (20) 

With the exception of the explicit three-body term Ai , this is the impulse approximation 
to the optical potential. 

Substituting (fTTj) into © gives 

t rr = Tr + Tr G bs BM r - (T R + T R G° BS BM R ) G° BS B (T PR + T 88 ) 
= T R + T R G° BS BM R - T R G° BS B (T PR + T RR ) 

-T R G BS BM R G BS B (T PR + T RR ) 
= T R + T R G° BS BM R - T R G° BS B (T PR + T RR ) + T R G° BS B (T RR - M R ) 
= T R -T R G BS BT PR . (21) 

This can be used in (JHJ) to yield 

T PR = -MpG° bs BT rr 

= -M P G BS B(T R -T R G BS BT PR ) 

= -M P G° BS BT R + MpG° bs BT r G° bs BT pr (22) 

Iteration of this shows that that it can be rewritten as 

t pr = _ T pp G ° bsBTr . (23 ) 

Using flU in (J2T|). 

t rr = Tr + Tr G° bs BT pp G bs BT r . (24) 
The complete t-matrix is the sum of (fT2"|) . (|2*3*|) and (|2"4"j) . This can be written as 

q- _ q-PP _j_ q-RP _|_ q-PR _|_ q-RR 

= 1 -l R Lr BS Lfl -1 Lr BS tSl R + 1 R + l R Lw BS til Lr B gJD 1 R 

= T R + (1 - T R G° BS B) T pp (1 - G° BS BT R ) . (25) 



III. WAVE FUNCTIONS 



We also need a similar separation for the scattering state of particle 1 with the bound 
state of particles 2 and 3 and for the three-body bound state. To obtain the former, consider 
the left-handed propagator defined by ([ A7|) . Using (|23|). this can be rewritten as 

Ql = G° BS - G° BS (1 + B) TG° BS 

= g bs ~ G bs (! + B ) t R G °bs 

-G° BS (1 + B) (1 - T R G° BS B) T pp (1 - G° BS BT R ) G° BS 

= G° BS - G° BS (1 + B) T R G° BS 

-G° BS (1 + B) (1 - T R G° BS B) iG-^D 1 \Y^) iQ^-iGoptKiQ 1 ) 
x(T^\D^(l-G BS BT R ) G BS 
= G° BS - G° BS (1 + B) T R G° BS 

+iG BS (1 + B) (1 - T R G° BS B) iG^D 1 \Y^) 

x (Gi - G^optGopt) iQ 1 <$ {2)1 | D lT (1 - BT R G BS ) . (26) 

From the residue of the pole contribution to G\, we can identify the scattering state as 

(S^l = {p u P^\{l-VoPTGoPT)iQ l {& 2)1 \D lT (1- BT R G BS ) 

= (p 1 ,P 1 \(l-T p T G 1 ) i Q 1 (^ 1 \D lT (1-BT R G° BS ) , (27) 

where P l is the momentum of the bound state of particles 2 and 3. 

The three-body bound-state vertex function, defined by (|A10J) . can be written as 

|r> = -VG° BS (1 + B) |r> = -MG BS B |r> . (28) 

Separating the vertex function into contributions where the last interaction contains either 
the pole in the residual parts of the scattering matrix for particles 2 and 3 gives 

|r P ) = -MpG° bs B |D (29) 
\Tr) = -M R G° BS B |D • (30) 

Using the definition of the pole contribution, 

|T P ) = D^T^iG^iQ^T^D^GlsBlT) 

= D^T^iQ'^^D^BlT) . (31) 



The remaining part of the vertex function is 

\r R ) = — (Tr + t r g° bs bm r ) g° bs b \t) 
= —TrG bs b \t) - TrG bs bja r g bs b |r) 
= —t r g° bs b |r) + t r g bs b \t r ) 

= -TrG bs B\T p ) (32) 

Using D 1T BD 1 = in (|3ip. substituting \Tr) from (|32[) and iterating once, 

|r P ) = D 1 \T® 1 )iQ x {$® 1 \D lT B\T R ) 

= —D 1 \T^) iQ 1 {^ x \D lT BT R G BS B \T P ) 

= —D 1 Ir^ 1 ) iQ 1 | D^BTrGIsBD 1 \Y^) iQ 1 | D lT B \T) 

= —D 1 Ir^iQ 1 {^ l \D lT BT R BD l (-zG^zQ 1 (S^ D lT B \T) 

= -1D 1 |r (2)1 ) V OPT (- iG i) < $(2)1 1 |r) • (33) 

Comparing the first and last lines shows that 

Q 1 ($ (2)1 1 D lT B \T) = -VoptHGx) ($ (2)1 1 D lT B \T) , (34) 

which is of the form of the bound-state vertex function for the optical model. 

The complete three-body vertex function can now be reconstructed using (J32*)) and (pUJ) 
to give 

|r> = |r P ) + \v R ) = (i - T R G BS B) \r P ) 

= (1-T r G bs B)D 1 \tW 1 )iQ 1 (^ 1 \D iT B\T) . (35) 
The Bethe-Salpeter wave function is then 

|*) = G° BS \T) = (1 - G° BS T R B) D 1 j^ 1 ) (-iG^iQ 1 {^\D lT B \T) (36) 

and 

^ OPT = (-iG^iQ 1 | D lT B |T) (37) 
can be identified as the optical model wave function. 
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(a) (b) 

FIG. 1: These Feynman diagrams represent contributions to the seven-point function. The particles 
are label 1 to 3 from top to bottom. The rectangles represent two-body kernels. 

IV. ELECTROMAGNETIC CURRENT MATRIX ELEMENT 

At this point it is necessary to deal with a problem that is occurs in describing the 
current matrix element for the Bethe-Salpeter equation that is not present in the usual 
nonrelativistic approach. Some care must be taken with defining the Bethe-Salpeter current 
operator if the matrix elements are to be of the form 

J»= (*/| J"|*i) . (38) 

Consider the contribution from the absorption of a virtual photon on particle 1. The initial 
state wave function can produce a contribution described by Fig. [TJt, while the final state 
wave function can produce a contribution described Fig. ^p. Since these are Feynman dia- 
grams and topologically equivalent, these two diagrams give identical identical contributions 
and including both contributions will double count. Therefore, in order to write the matrix 
element element in the symmetric form (J38j) . it is necessary to correct the current operator 
to correct for the double counting. This can be done by replacing the one-body current by 



40] 



39, 401 and is included in the 



the currents represented by Fig. |21 This was pointed out in 
definition of the Bethe-Salpeter effective current operator defined in |4l| given by (jA18|) and 
i"23j). 

Now consider the electromagnetic current matrix element for ejecting particle 1 from the 
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FIG. 2: Feynman diagrams representing the correction to the current operator to correct for double 
counting. 



bound state into the continuum state where particles 2 and 3 remain bound. This is 
= (p^P 1 ! (1 -Topri-iG^iQ 1 (^\ D lT (l - BT R G° BS ) 



where 



(1 - G° BS T R B) D 1 \^) {-iG^iQ 1 {^ l \D lT B \T) 
= (p lj P 1 \(l-T PTHG l ))J£ PT (-iG 1 )(^ 1 \D lT B\r) , (39) 

r oPT = ^ <$ (2)1 | D lT (1 - BT R G° BS ) J? s (1 - G° BS T R B) D 1 feW^iQ 1 (40) 

and is defined by (|A23J) . 

Considerable care must be taken in evaluating this expression. To simplify the derivation, 
we have used the operator iQ} to place the bound state on shell. Operators of this type 



were introduced in [ 381 and elaborated in 



3 and H. 



This is a very singular operator 



and must be treated with extreme care. In particular this operator is not associative and 
its evaluation depends on its context in the evaluation of physical quantities. To see this 
consider the two-body scattering matrix for particles 2 and 3 given by 

M 1 = V 1 — V l G l M l = V l - M x G l V l . (41) 

The second form of this equation can be solved to give 

V 1 = M 1 + M l G 1 V 1 . (42) 

12 



Substituting this into the first form leads to the nonlinear form of the equation for the 
scattering matrix 

M 1 = V 1 - M 1 G 1 M 1 - M 1 G 1 V 1 G 1 M 1 = V 1 - M 1 G 1 G 1 ~ 1 G 1 M 1 - M 1 G 1 V 1 G 1 M 1 



V 1 -M l G 1 {G 1 1 + V 1 )G 1 M 1 = V l - M 1 G 1 Q 1 l G 1 M 



(43) 



Note that since Q 1 1 must vanish at the bound state pole, both sides of this equation have 
simple pole at this point. The residues of these poles give 

-jr^^zQ^r^ 1 ! =-|r( 2 ) 1 )2Q 1 <r( 2 ) 1 |G 1 ^ 1 " 1 G 1 |r( 2 ) 1 )zQ 1 <r( 2 ) 1 | . (44) 

This requires that 

iQ 1 = iQ l (r (2)1 | G 1 Q 1 ~ 1 G 1 ll^ 1 ) iQ 1 = iQ 1 G 1 ' 1 \^) 1Q 1 . (45) 

Clearly, if we choose to group Q 1 1 with either the first or last occurrence of Q 1 on the right 
had side of this equation, the right-hand side will vanish and the equation will be violated. 
This means that the operators on the right-hand side must be evaluate as a whole without 
attempting to evaluate them in a pair- wise manner when appearing in this context. 

As a first step in simplifying the optical model current operator, consider the effective 
current define in (IA23I) as 



'eft' 



j(i)m (1 + B) + (1 + B) J{L (1 + B) = J^DD T + DD T J»DD T 



= D{J^ + J^ + J^ + J^ + J^ + J^ + J^D T 

= D ( J 1 ^ + iV 1 Jf + J 2 » + + + J% + + J*£) D 1 

where we have used the identity 1 + B = DD T with 

/l\ 
1 



(46) 



D 



1 



(47) 



Consider 

J 1 ^ + Jj£ = + U?V l + J\ x (48) 

where J\ x is the two-body current for particles 2 and 3 and the remaining terms contain 
only one-body currents. The one-body currents can be rewritten as 

J 1 " + iV 1 Jf = U^G 1 ' 1 + iV x J{ = iJ^g 1 ' 1 , (49) 
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Since these pieces contain the inverse of the interacting propagator for particles 2 and 3, 
care must be taken when these contributions are associated with the operator iQ 1 . For this 
reason we need deal with the contributions of these currents to the optical model current 
separately. This gives 



i _ 

= i 
-i 



iQ 1 (^j D lT (1 - BT R G° BS ) DiJ^Q 1 - 1 D T (l - G° BS T R B) D 1 \^) tQ 1 

Q 1 (^\ iJ^g 1 ' 1 \^ 1 )iQ 1 -iQ 1 {^ l \iJ^g 1 - 1 D T G BS T R BD l \^ x )iQ l 
Qi | D lT BTRGzsDU^g 1 - 1 \<$>^ l )iQ l 
+tQ^^ 1 \D lT BT R G° BS DzJ^g 1 - 1 D T G° BS T R BD 1 \^ 1 )tQ 1 
= iJ^ + iQ 1 <$ {2)1 | D lT BT R G% s DiJ^g 1 ' 1 D T G% s T R BD l \^)iQ l . (50) 

This has been simplified using f)13j) and the identities 

iQ l g l ^O = Og l '\Q} = 0, (51) 

where O is any operator other than iQ} or g l . The complete optical model current 
operator therefore reduces to 

r QpT = UfiQi + iQ} <$( 2 )!| [D lT BT R G BS D i J^g 1 - 1 D T G BS T R BD 1 
+D lT (1 - BT R G° BS ) D ( J 2 " + J 3 " + Jffi + J% + + Jg) 
xD^l-GMD 1 ]!^ 1 )^ 1 . (52) 

Note that the first term yields the usual impulse approximation while the remaining contri- 
bution corresponds to a considerable number of Feynman diagrams. 

To show that the optical model current is conserved we need the identities 

g„Jf = [ ei (q),G^] (53) 

<bJ% = [e2(q)+e 3 (q),V 1 ] (54) 

%J iil = [ei(g),G^s _1 ] (55) 

<bJ& = [er(q),V] (56) 

where ei{q) combines the charge operator for particle i with a four- momentum shift operator 
and 

er(q) = e^q) + e 2 (q) + e 3 (q) . (57) 
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Using these identities, contraction of the optical model current with the four-momentum 
transfer gives 

%Jopt = i[e x {q),G^]iQ l + [D lT BT R G° BS Dt[e 1 (q),G^}Q 1 - 1 D T G BS T R BD 1 
+D lT (1 - BT R G° BS ) D([e 2 (q) + e 3 (q), G^ 1 } + [e 2 (q) + e 3 (q), V 1 ] 
+ [e T (q), V° + V 2 + V 3 ])D T (l - G° BS T R B) D 1 } \<S>^) iQ l 
= ileM.G^iQ^iQ 1 <$ {2)1 | [D lT BT R G BS Di[e l {q) 1 G^]g 1 - 1 D T G BS T R BD 1 
+D lT (1 - BT R G° BS ) D([e T (q), G ^' 1 ) + [e T (q), V° + V 1 + V 2 + V 3 ]) 
xD T (l-G° BS T R B) D 1 

— D lT (1 - BT R G° BS ) D[ ei (q), G^' 1 + V']D T 

x (1 - G° BS T R B) D 1 } fiV^iQ 1 (58) 

Using 

zG^g'-^G^ + V 1 , (59) 

this can be simplified as 

%Jopt = 2[ei(g),G^W+*2 1 ($ (2)1 | [D lT (1-BT R G BS )D 

x (ie T (q), G^ 1 ] + [e T (q), V° + V 1 + V 2 + V 3 ])£> T (l - G° BS T R B) D 1 
-[e^^ls- 1 + V 1 ] + [e^G^T 1 + V^D^^BD 1 
+D lT BT R G BS D[e 1 (q),G BS ~ 1 + V 1 ]] l^zQ 1 . (60) 

From (|5I|I . 

iQ^ei^^^+VW^O. (61) 

So, 

g M J^ pT = i^qlG^liQ' + iQ 1 ^ 2 ^ [D lT (1 - BT R G° BS ) D 

x ([e T (g), G ^ 1 ] + [e T (q), V° + V 1 + V 2 + V 3 ])£> T (l - G° BS T R B) D 1 
+ei(q) (cPes' 1 + V 1 ) D T G° BS T R B D l 

-D lT BT R G° BS D fag- 1 + V 1 ) ei (g)] ^W^iQ 1 . (62) 

No consider 

(GV 1 + V 1 ) D T G° BS T R BD 1 = D lT [(l + S)T Rj B + V(l + B)G° BS T R B]D 1 . (63) 
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This can be simplified using 

VG° BS T R B = VG° BS M R (1 - G° BS BT R ) B 

= VG° BS (M - M P ) (1 - G° BS BT R ) 
= (V - M + M P ) (1 - G° BS BT R ) B 
= (V - M R ) (1 - G° BS BT R ) B 

= VB - VG° BS BT R B - T R B (64) 

such that 

(g bs ~ 1 + V 1 ) D T G° BS T R BD 1 = D lT [(l + B)T R B + VB - VG° BS BT R B - T R B 

+VBG° BS T R B]D 1 
= D lT [BT R B + VB]D 1 = D^BTrBD 1 . (65) 

Similarly, 

D lT BT R G° BS D (els' 1 + V 1 ) = D lT BT R BD l . (66) 

Note also that 

D[e T (q),G BS ' 1 + V° + V 1 + V 2 + V 3 }D T = e T (q)(l + B)^' 1 + V(l + B)] 

-\G BS - 1 + (1 + B)V](1 + B)e T {q) , (67) 

and 

[Gps' 1 + V(l + B)] (1 - G° BS T R B) D 1 

= [G BS ^ + V(l + B) — T R B - V(l + B)G° BS T R B]D 1 

= [G ^ 1 + V(l + B) — T R B -VB + VG° BS BT R B + T R B - VB) G° BS T R B] D 1 
= [G^- 1 + V}& = Duels' 1 + V 1 ] . (68) 

Similarly, 

D lT (1 - BT R G° BS ) [G BS - 1 + (1 + B)V] = [G° 5 _1 + V l ]D lT . (69) 
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Using these identities we can rewrite 

+tQ} (<S>^\ [D lT (1 - BT R G° BS ) e T (q)(l + B)D 1 [G° BS ~ 1 + V 1 ] 
-[Gls' 1 + V 1 ] J D lT (l + B)e T {q) (l - G° BS T R B) D 1 
+e 1 (q)D lT BT R BD 1 - D lJ BT R B D 1 ei (q)] |$ (2)1 > tQ 1 

= i[ ei (q), G^iQ 1 + e^q^Q 1 (& 2)1 \ D lT BT R B D 1 \^) tQ 1 
-tQ 1 (^\ D lT BT R BD l l^ 2 ) 1 ) iQ 1 e 1 (q) 

= ileMiG^iQ 1 + [eM^Vopr] (70) 

Or, 

% j opt = A^M^opt] ■ ( 71 ) 

This is the Ward-Takahashi identity for the optical model current and along with the wave 
equations for the optical model wave functions guaranties that the current matrix elements 
will be conserved. 

V. CONCLUSIONS 

We have shown that the three-body Bethe-Salpeter equation can be reduced to an effective 
two-body optical model. An effective current appropriate to this model has been constructed. 
This current is shown to satisfy a Ward-Takahashi identity involving the optical potential 
which results in conserved current matrix elements. This conserved current contains a 
substantial number of contributions not included in current RDWIA calculations and the 
contributions of these extra terms in various kinematical regions need to be considered 
carefully. 

Although this paper deals with a simple three-body system, extension of this approach 
two include additional constituents is possible as will be described in a subsequent paper. It 
may also be useful to consider limiting cases of this approach to understand its relationship 
to the mean field approaches used for most calculations. 
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APPENDIX A: REVIEW OF THE MATRIX FORM OF THE BS EQUATION 
FOR DISTINGUISHABLE PARTICLES 



This appendix contains a short summary of the matrix form of the three-body Bethe- 
Salpeter equations and effective current as defined in Q. 

The three-body Bethe-Salpeter equation can be obtained by examining the sum of all 
Feynman diagrams contributing to the three-body scattering matrix. Contributions to these 
diagrams can be classified according to whether the contribution can be separated by cutting 
only the three propagators associated with the external legs of the scattering matrix. Those 
diagrams which can not be separated in this way are three-body irreducible diagrams. The 
irreducible diagrams fall into two classes: those where only two of the three particles are 
interacting and those where all three particles are interacting. The sum of all three-body 
irreducible diagrams is represented by the kernel V° and the two-body irreducible diagrams 
contribute to the two-body kernels V 1 where only particles j and k (with j ^ k ^ i) are 
interacting. 

The complete scattering amplitude can then be written in terms of an integral equation 
with the above mentioned kernels. It is convenient to express the complete scattering matrix 
in terms of subamplitudes where the indices i and j characterize the subamplitudes 
according to the character of the last and first interactions; that is, for i ^ (j 7^ 0), the 



particles i (j) are not taking part in the last (first) interaction, i — (j — 0) means 
last (first) interaction is genuine three-particle interaction, equations (see also ref. 

3 



that the 
): 



T ij = V%j - VG° BS J2 Tkj , ( Al ) 

k=0 

where 

( V° for i = 
V 1 = { (A2) 
[ VHGJ 1 for i = 1,2,3 

and G° BS = — GiG 2 G 3 . The form of these equations suggests that it is convenient to represent 
this set of equations in a matrix form. Defining the matrices (V)«j = V l 5ij, (T)ij = T l i and 
(B)ik = 1 — 5ik for i, j = 0, 1, 2, 3, the three-body scattering equations can be written as 



T = V - VG BS (1 + B)T = V - T(l + B)G BS V , (A3) 



where G° BS = G° BS 1. 
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Numerical solution of these integral equations requires that they must be put in a form 
where the kernels are connected or can be made to be connect by iteration. This is done by 
reexpressing the equations in terms of two- and three-body t-matrices defined in terms of 
the corresponding interaction kernels as: 

M = V - VG° BS M = V - MG° BS V . (A4) 

The complete t-matrix can then be written as 

T = M - MG° BS BT = M- TBG° BS M . (A5) 

In matrix form, it is necessary to define right- and left-handed propagators 

G R = G BS - G° BS T (1 + B) G° BS = G° BS - G° BS V (1 + B) Q R (A6) 

and 

G L = G BS - G° BS (1 + B) TG° BS = G° BS -Q L (1 + B) VG° BS . (A7) 
The inverses of these propagators are 

Q-r = (G BS )- l + V(l + B) , (A8) 
Qt = (G° B5 )- 1 + (l + B)V. (A9) 

The bound state can be obtained from consideration of the singularities of the t-matrix. 
and satisfies the equations 

|r> = -VG BS (1 + B)\T) , (A10) 

<r| = -(r|(i + B)G° B5 v. (ah) 

Defining the Bethe-Salpeter wave function as |\I>) = G° BS \T) these can be rewritten as 

[(CPbs)- 1 + V (1 + B)} |*) = g R l I*) = , (A12) 
(*l [(GbsY 1 + (1 + B)V] = <*| Gl 1 = . (A13) 

The scattering wave functions also satisfy the same equations. 

The three-body Bethe-Salpeter current can be determined by considering all diagrams 
contributing to the seven-point function with six legs corresponding to the three incoming 
and outgoing particles and one photon leg. By separating the diagrams into three-particle 
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reducible and irreducible contributions the current operator can be identified as the sum of 
all irreducible seven-point functions. There will be three types of contributions: one-body 
contributions where the photon attaches to only one of the interacting particles, two-body 
contributions where the photon attaches internally to a two-body interaction, and three- 
body contributions where the photon attaches internally to a three-body interaction. 
The one-body currents are of the form 

= -JfG^Gj 1 , (A14) 

where t = 1, 2, 3, r ^ s ^ t and Jj* is the vertex for attaching a photon to particle t, for 
which q^Jf = [et(q),G^ ]. Each of these currents satisfies the Ward-Takahashi identity 

^= [etiqKGBs)- 1 } • (A15) 

The two- and three-body currents follow from attaching a photon line to all particle lines 
and into momentum-dependent vertices internal to the two- and three-body Bethe-Salpeter 
kernels. Defining Jrs^ as the two-body current associated with the two-body contributions 
to the kernel and j( 3 )^ as the completely connected three-body current, one can write the 
exchange current for the three-body system as 

\ J^iGT 1 , t = 1,2,3 (r^s^t) 
J%={ n (A16) 

{ J(3)M t = 

which satisfies the Ward-Takahashi identities 

T tu J Mg) + e s (g),V], t = 1,2,3 (r ^ s ? t) 

^ J ex = < r m (A17) 

{[e T (q),V°], t = 

where e T (q) = e x {q) + e 2 (g) + e 3 (q). 

Following the argument in Section IIV| it is necessary to include a contribution to the 
effective current to compensate for the double counting inherent in the symmetric expression 
for the current matrix element. This can be done by defining the interaction current 

Jfi = J% + iV f J? for t ? , (A18) 

Jfi = >% = J {3) " ■ (A19) 

The matrix form of the effective current is obtained by first defining the total one-body 
current as 

3 

jW" = J* , (A20) 



i=i 
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and denning a diagonal matrix with components denned by (|A18|) and (|A19[) : 



int 



(Hag (J in t, J^f> ^int' ^mt) 



% J L = M<?),v] . 



The effective current can then be identified as 



jm = j(Dm {1 + B) + {1 + B) j L (1 + B) 



Contraction of the four-momentum transfer with the effective current gives 

erOz), Cls' 1 ] (1 + B) + (1 + B) Mg), V] (1 + B) 
e T (g)(l + J B)a i? 1 -a L 1 (l + B)e T (g). 



<?M J cff 



(A21) 
(A22) 



(A23) 



(A24) 



So, using the wave equations, the current will be conserved. 
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